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Abstract 

It is well known that the two-dimensional boundary layer on a concave wall 
is centrifugally unstable with respect to vortices aligned with the basic flow 
for sufficiently high values of the Gortler number. However, in most 
situations of practical interest the basic flow is three-dimensional and 
previous theoretical investigations do not apply. In this paper the linear 
stability of the flow over an infinitely long swept wall of variable curvature 
is considered. If there is no pressure gradient in the boundary layer it is 
shown that the instability problem can always be related to an equivalent two- 
dimensional calculation. However, in general, this is not the case and even 
for small values of the crossflow velocity field dramatic differences between 
the two and three-dimensional problems emerge. In particular, it is shown 
that when the relative size of the crossflow and chordwise flow is 0(R ), 

where R is the Reynolds number of the flow, the most unstable mode is time- 
dependent. When the size of the crossflow is further increased, the vortices 
in the neutral location have their axes locally perpendicular to the vortex 
lines of the basic flow. In this regime the eigenfunctions associated with 
the instability become essentially 'centre modes' of the Orr-Sommerfeld 
equation destabilized by centrifugal effects. The critical Gortler number for 
such modes can be predicted by a large wavenumber asymptotic analysis; the 
results suggest that for order unity values of the ratio of the crossflow and 
chordwise velocity fields, the Gortler instability mechanism is almost 
certainly not operational. 
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1. INTRODUCTION 


In recent years there has been much interest in the effect of boundary 
layer growth on the instability mechanisms to which boundary layers are 
susceptible. Our concern here is with the Gortler vortex instability 
mechanism which is known to occur in two-dimensional boundary layer flows over 
concave walls. The latter problem has been discussed by Gortler (1940), 
Hammerlin (1955, 1956), Smith (1955), Floryan and Saric (1979) and Hall 

(1982a, b, 1983). In the latter three papers it was shown that the parallel 
flow approximation used previously to reduce the linear stability equations to 
ordinary differential equations is justifiable only for vortices of small 
wavelength. In this case the effect of boundary layer growth can be taken 
care of in a self-consistent manner and an asymptotic form for the right-hand 
branch of the neutral curve can be derived. At order one values of the vortex 
wavelength the parallel flow approximation is not valid and the full 

linearized partial differential equations must be solved numerically. The 
reader is referred to the paper by Hall (1983) for a discussion of the results 
and implications of such a calculation, it suffices here to say that the 
physically unacceptable results of previous calculations in this wavenumber 
regime are shown to be a direct consequence of the parallel flow 

approximation. 

There are many practical situations where Gortler vortices are thought to 
occur but in most of these the basic boundary layer flow is three- 

dimensional. If, for example, we are concerned with the flow on the concave 
region of a swept laminar flow wing or that over a turbine blade, then the 
three-dimensional nature of the basic flow is important and cannot be 
neglected; thus previous theoretical stability calculations do not apply. In 
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order to determine the effect of three-dimensionality on Gortler vortices for 
such flows we shall in this paper investigate the instability of the boundary 
layer on the concave part of an infinitely long swept cylinder. 

The basic boundary layer flow associated with this flow is easily 
calculated using the method of Sears (1948). The parameters which 
characterize the flow are a, the angle of yaw, R the Reynolds number and 6 
the curvature parameter. If ot = 0, we know that steady Gortler vortices 

occur for 0(1) values of R^ 2 6 and are aligned with the streamlines of the 
basic flow. The aim of this paper is to determine how the structure of the 
instability changes when a is increased from zero. The present analysis is 
restricted to small values of the vortex wavelength, but our experience with 
the two-dimensional problem would suggest that the results which we obtain are 
useful at 0(1) values of the wavelength. However, the differences between 
small wavelength Gortler vortices in two- and three-dimensional boundary 
layers are so fundamental that we expect that vortices with 0(1) wavelength 

cannot even occur in a three-dimensional boundary. 

The relative size of the crossflow and chordwise flow is the crucial 
factor in determining the structure of Gortler vortices in three-dimensional 
boundary layers. Suppose that X is a typical value for the ratio of the 

latter velocities in the boundary layer. We show that, when X is increased 

from zero, the first crucial change in the structure of the instability occurs 
when X ~ R~ ^ 2 . In this regime the vortices become time-dependent and 
meander as they develop in the chordwise direction. Surprisingly , the 
orientation of the most dangerous mode is always determined by the vortex 
lines of the basic flow. We show that at the neutral location the vortex 

boundaries align themselves so as to be locally perpendicular to the vortex 
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lines of the basic flow. Thus, in general, unlike the two-dimensional case, 
the vortices are not locally parallel to the flow direction. At larger values 
of X further dramatic changes in the vortex structure emerge, until finally 
the eigenfunctions develop into 'centre modes' of the Orr-Sommerfeld equation 
made unstable by centrifugal effects. We are able to find the development of 
the vortices up to the regime X ~ R ^ and it appears that for X » R 
the Gortler vortex instability mechanism is either not present or cannot be 
described by asymptotic means. Thus, we are able to describe the evolution 
of Gortler vortices for 0 < X < 0(R -1 / 8 ) but it turns out that by writing 
X = Xr - 1 ^ 2 all the regimes of interest can be obtained by first letting 
R with X fixed and then taking further limits on X. This is the 
procedure which we will adopt in this paper. 

We find that the above results do not apply when the basic three- 
dimensional boundary layer has zero pressure gradient* In this case, the 
stability equations can be solved for 0(1) values of ot and, in fact, can 
be reduced to an equivalent two-dimensional problem. The vortices are then 
aligned with the streamlines of the basic flow again and the most dangerous 
modes are steady. 

The procedure adopted in the rest of this paper is as follows. In Section 
2 we formulate the stability equations for the flow over a slightly yawed 
infinite cylinder and solve these equations asymptotically for small vortex 
wavelengths. In Section 3 we discuss the special case of flows with zero 
pressure gradient and in Section 4 we discuss our results and their practical 
implications. 
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2. GORTLER VORTICES IN A SLIGHTLY THREE-DIMENSIONAL BOUNDARY LAYER 

We consider the stability of the boundary layer flow over the cylinder 
y = 0, -<» < z < <». The z axis is taken to be a generator of the cylinder 
and y measures distance normal to the surface* The x coordinate measures 
distance along the curved surface which is taken to have variable curvature 

| X 

— < (j) where a and Z are length scales. The Reynolds number R, the 
curvature parameter 6 and the Gortler number G are defined by 

V 

R = -T r » (2.1a) 


6 



(2.1b) 


G 



(2.1c) 


where U Q is a typical velocity in the x direction. We assume that the 
Reynolds number is large, whilst 6 is taken to be small. More precisely, we 
consider the limit 6 -*■ 0, with G held fixed. The basic three-dimensional 
boundary layer flow is assumed to be of the form 


u = U 0 (u(X,Y),R" 1/2 v(X,Y), Aw(X,Y))(1 + 0(R 1/2 ), 


where 




1/2 

> 


and A is to be specified shortly. Following the procedure of Sears (1948) 
u and v are found by integrating numerically the two-dimensional boundary 
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layer equations. Xw is then calculated from the z momentum equation. At 
this stage we need only insist that the basic flow be independent of the 
spanwise coordinate, later we shall indicate the relevance of our calculation 
to some particular boundary layer profiles of practical importance. 

The appropriate scaling for X, which, fixes the angle of yaw of the 
cylinder, can be found by taking X ~ , y > 0 and varying y until 

the Gortler instability mechanism driven by u is modified by the crossflow 
at zeroth-order. It is known from the two-dimensional problem that the 

characteristic wavelength of the vortices is of the same order of magnitude as 

g 

the boundary layer thickness, so that the convective terms Uq u ^ and 

3 1/2 

U 0 Xw will be of comparable order when R * 00 if X ~ R • We write 


X 



to effectively restrict our attention to flows over cylinders yawed at an 
- 1/2 

angle a ~ 0(R ) to the oncoming flow* At a later stage we can consider 

the limits X + 0, X + ® in order to recover some information about the 
- 1/2 - 1/2 

regimes X « R and X » R respectively* We define the variable Z 

by 


Z 


R 1/2 


z/£ 


and perturb the basic flow by writing 

u = U Q (u + U(t,X,Y)E, vR -1/2 + V(t,X,Y)R _1/2 E, Xr” 1/2 w 

+ W(t,X,Y)R -1 ^ 2 E)(l + 0(R -1 ^ 2 )) , 


(2.1c) 
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where t is a time variable scaled on £/Uq whilst 

E = exp( iaZ) . 

We then substitute the above expression for u_ into the Navier-Stokes 
equations and linearize to obtain 

U t + aU X + Ua X + ^ U Y + v “y + ^ ”" iaU = U YY “ a2 U> (2.2a) 

V + uV^ + Uv^ + vV^ +Vv^ + X wiaV + Gk(X)uU = -P^ + V> (2.2b) 

+ uW^ + X Uw x + vW^ + X Vw Y + X wiaW = -iaP + W YY W, (2.2c) 

U v + V v + iaW = 0. (2. 2d) 

Here P is the nondimensional pressure perturbation corresponding to the 
disturbed velocity field and we have neglected terms of relative order 
The equations (2.2) reduce to the corresponding two-dimensional equations when 
X =0. We shall see later that if u = w we can solve the three-dimensional 
stability problem even for X = 0(1). In fact, in this case the three- 
dimensional problem can always be reduced to an equivalent two-dimensional 
problem so that the results of Hall (1982a, 1983), hereafter referred to as I, 
II, can be used. For the remainder of this section we assume that IT t w and 
determine how the crossflow modifies the instability. 

It is clear that if the wavenumber a is 0(1) then (2.2) must be solved 
numerically in the manner described in II, but we must allow for a possible 
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time dependence of the instability. Our experience in I, II showed that it is 
only in the limit a * 00 that the usual results of stability theory apply 
and, in fact, for a~ 0(1) there is no such thing as a unique neutral 

curve. However, we found that the asymptotic results for a >> 1 gave useful 

information about the instability for a~ 0(1), so we shall concentrate here 
on the further limit a -*■ 00 . 

It was shown in I that for a >> 1, the vortices concentrate themselves in 

an internal viscous layer of thickness In the neutral location this 

layer is in the position where u(X,Y) most violates Rayleigh's criterion. 
If the vortices are not locally neutrally stable, this is not the case and we 
shall show that for sufficiently large crossflow velocities the position of 
the viscous layer where the vortices are concentrated is never fixed by 

Rayleigh's criterion* 

We are now in a position to determine the asymptotic solution of (2.2) for 
large values of the nondiraensional wavenumber a. Before doing so, we must 
specify the relative size of X, which determines the magnitude of the 
crossflow velocity field. The aim of our calculation is to determine the 

effect of increasing X on the Gortler vortex instability mechanism, so it is 
convenient for us to take the limit a + <», X~ a J for J > 0. Afterwards 

we shall indicate for what regimes of the curvature and flow velocities these 

further limits apply. It turns out that there are three major regimes for 
J which enable us to describe how a Gortler vortex eigenfunction evolves into 
essentially an Orr-Sommerf eld eigenfunction when the crossflow is increased. 
These three regimes correspond to (a) a~ X , (b) a~ X ' , and (c) 

j /3 

a ~ X , respectively. We shall see that for T » 1 , the first significant 

9 

difference between the two- and three-dimensional problem emerges when a ~ a • 
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If a is decreased to 0(X” 3 ^ 3 ) the mechanism which enables the vortices to 
be concentrated in an internal viscous layer changes. Thus for the two- 

dimensional problem the decay of the vortices is due to the fact that the 
local Gortler number has a maximum at the centre of the layer whereas for the 

o /c 

three-dimensional problem with a ~ X , the decay of the vortices is 

—1/3 

facilitated by convective effects. A further decrease in a to 0(X ) 

shows that the Gortler vortices at this stage are essentially Orr-Sommerfeld 

1/3 

eigenfunctions destabilized by centrifugal effect. Moreover, for a a ~ X 

a neutral curve with both left and right hand branches and a minimum Gortler 

_ 1/3 

number is obtained, this suggests that for a a << X , instability is almost 
certainly impossible at finite values of G. We now describe the three 

regimes described above: 

a. The limit a + 00 with X ~ 0(a)^ 2 , G ~ O(a^). 

It is known from the two-dimensional problem that for a » 1 the flow is 

neutrally stable for G ~ a^. We anticipate that this will also be the case 
— 1/2 

for X ~ a and therefore write 

G = g Q a 4 + a 3 + •••, (2.3) 

— 3 * 1/2 

whilst X and ■%— are replaced by X a and ia with 

O t 

a - a 3 ^ °o + aa l + **** (2.4) 

The aim of our calculation is to determine the frequency o and Gortler 
number G such that the flow is locally neutrally stable at X. However, in 



9 


order to show how non-neutral disturbances evolve with X, we allow for 
growth or decay in X by writing 

U = I a n/2 U n (X,n)exp{iot + a 2 f 0 Q (X) + a‘ 1/2 ^(X) +,*-dX} , (2.5) 

together with similar expansions for V/a 2 , W/a 3 / 2 , P/a 5 / 2 . The variable n 
is defined by 

n = {y - Y(X)}a 1/2 , (2.6) 

where Y = Y(X) is the location of the viscous layer where the vortices are 
concentrated. The thickness of this layer was shown in I to be 0(a -1 / 2 ) 
and it corresponds to a second-order turning point of a WKB solution of 
(2.2). The growth rate functions 3 (X), i=0,*«*, appearing in (2.5) are in 
general, complex quantities and for given G and o will clearly vary with 

X. The velocity components of the basic flow expand locally around Y = T(X) 

as 

u = U 0^ X ) + a Uj(X)q + a * u 2 (X)n 2 +,*•*, 

together with similar expansions for v and w. It is a routine matter to 

substitute the above expansion into (2.2) and by successively equating like 
powers of a “l/ 2 we generate a sequence of equations to determine 
( U n ,V n’ W n’ P J f ° r n = ^ system corresponding to n = 0 is 

(u 0 3 0 + 1)U Q = -V Q uj, (2.7a) 


^ u 0 ^0 + ^0 *’ S 0 KU 0 U 0 ’ 


(2.7b) 
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w 


0 



(2.7c) 


p o * 0 + u o e ol lw o- 


(2.7d) 


At this order the crossflow has no effect on the expansion procedure and 
(2.7a,b) have a solution if 


( u 0 e o + ^ g o u o u 


1 


( 2 . 8 ) 


so that instability occurs if u^ u^ k > 1 and the growth rate 8 q is 

real. The functions Vq, Wq and P Q can then be expressed in terms of Uq 
using (2.7a,b,c). At next order the system of equations which determines 
(Ul V}, W^, P^) is found to be 


[«0 6 0 + *]»1 +»! V, --n[u, 6 0 u 0 + 2 u 2 V 0 ] 


- iX “o U 0 ■ io 0 U 0 " u 0 6 1 “o’ 


[»0 B o + llv i + g 0 Ku 0 “l ° ■ ”K e o v o + s o KU 1 “ol 


- iXw 0 V 0 - io 0 V 0 - u 0 6j V 0 , 


W 1 - 1 V l„ - i6 0 “o’ 


Pj - iWj[l + u 0 6 0 ] - 1W 0 Uj 8 0 0 - X» 0 H 0 - 0 Q W 0 + iu 0 0! V 


The above system has a consistent solution if 
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and 


ia Q + u Q + iXw Q = 0, (2.8a) 

23 o ujl + B 0 u Q ] = g Q k[2u q u 2 + u 2 ] . (2.8b) 


The first of these equations shows that is purely imaginary and for a 

given value of a ^ is determined by 


3 


li 



(2.8c) 


so that the orientation of the vortices varies with X. The second of the 
above constraints fixes Y = Y(x) as the location of the layer. In fact if 
gQ is given we see from (2.7), (2.8b) that at any value of X the growth rate 
and Y are determined by solving 


F(X,Y) = 0, F y (X,Y) =0, Y = Y(X) , 

where 

F = (u8 q + l) 2 - g Q u u y k. 

Thus, at any value of X the position of the viscous layer adjusts to produce 
a local maximum (in Y) for 8 q. In the neutral case, 8 q = 0 and the 
position of the viscous layer corresponds to (u u^) = 0, which is, of 

course, where Rayleigh's criterion is most violated. We note that at this 
order [Uq, Vq, Wq, Pq] and [Uj, Vj, , P^] remain undetermined, so that 

it is necessary to proceed further. In fact the next order system of 

equations can be solved if a consistency condition is satisfied, and this 
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leads to the required differential equation for Uq. We thus obtain 


3 2 U Q 2n{3 x n l + Wj Xi}U Q 
3n 2 {3 + u 0 e 0 } 


{g 2 u Q uj < - 2(1 + u Q 3 0 )(ia 1 + u Q 3 2 )}u Q 


{! + u 0 3 0 H3 + u Q 3 0 1 


- hn U Q = 0, 


(2.9) 


where 


h = 


f yy (x,y) 


2(1 + 


"o 6 0 ,(3 


u o s o* 


We are interested in the solutions of (2.9) which decay when n and 

after some manipulation we find that the most rapidly growing solution (in 
the X direction) can be written 


U 



{n + 


i[e u Uj + Wj X] 2 

h l 3 + u 0 3 0 J 1 


( 2 . 10 ) 


Here we have used the fact that & lr is zero and the eigenrelation 
corresponding to this eigensolution is 


{g 2 Vi < - 2[i + u 0 e 0 ]u 0 e 2r } 

{1 + u 0 + u o 3 o } 


[e u Ul + w x x ] 2 _ /r . 

o i 

h[3 + u Q 0 Q ]^ 


( 2 . 11 ) 


°1 + u 0 3 2i 


0. 


These equations determine 3 2r , 3 2i as ^ unct ^ ons °f X as the disturbance 
develops in the X direction. If we are interested in the neutral case, we 
can set 3 q = 3 2r =0 in (2.7), (2.11) to obtain the following two-term 
asymptotic expansion of the Gortler number, 
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1 L o U, + w. X} 

G = - {a 4 + 3a 3 [/hT + — — — — - ] + ••• } , (2.12a) 

where hg represents h evaluated with 8 q = 0 whilst 3^ is determined 
by (2.8c). The three-dimensionality of the basic flow affects the Gortler 
number through the squared term in the curly brackets and therefore has a 
stabilizing influence on the basic flow with X = 0. 

However, if the frequency <Jq is chosen such that 

U 1 ^ li + * = 0, (2.12b) 

.. O 

the Gortler number, correct to 0(a), reduces to its value for the flow with 
X = 0* Thus, even at this early stage, the importance of the modes which 
satisfy (2.12b) when they are locally neutrally stable is apparent. The 
constraint (2.12b) requires that the boundaries between neutral Gortler 
vortices be locally orthogonal to the vortex lines of the basic flow. 

We further note that steady Gortler vortices (a Q = 0, 6^ * 0) and 
Gortler vortices propagating in the spanwise direction (o^ t 0, 0^ = 0) are 
also possible. The former modes are locally parallel to the streamlines of 
the basic flow and are more stable than the oblique modes satisfying (2.12b). 
We shall see in the remainder of this section that it is the modes 
corresponding to (2.12b) which are important for X* » a , indeed, for X 
sufficiently large, the first-order term in (2.12a) can even be modified by 
three-dimensional effects. 
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— 5/3 4 

b. The limit a ->• 00 with X ~ 0(a ), G ~ 0(a ) 

1/0 

We have seen above that for X ~ 0(a ' ) the decay of the vortices away 

from the centre of the viscous layer is facilitated by the local maximum of 

the effective Gortler number at that level. We shall now show how this decay 

is facilitated by convective effects when X is increased further. In order 

to simplify the details of our expansion procedure we shall concentrate our 

attention on determining the neutral locations of the most dangerous Gortler 

— 1/2 

vortex modes. We have seen above that even for X ~ a the most dangerous 
modes are locally perpendicular to the vortex lines of the basic flow, and it 
is the development of these modes which we consider below. 

Suppose then that we concern ourselves with the mode which is locally 
neutrally stable at X = X n with a viscous layer located at Y = Y(X n >. The 

crucial property of the most dangerous three-dimensional modes discussed in 

0 0 

(a) is that at the viscous layer they cause the operator ^ + u + X aiw 

2 — 

to expand locally as sXi(Y — Y) where s is a constant. If X is taken to 

2 

be 0(a), this causes a term proportional to it) to appear in the equation 
corresponding to (2.9). This term, which is of course due to convective 
effects, enhances the decay of the vortex away from the viscous layer. In 

* C / O 

fact, if X = Xa , the decay of the vortex away from the viscous layer is 
entirely due to convective effects. However, at this stage the expansion 
procedure of (a) must be significantly altered because the thickness of the 
viscous layer now decreases to O^ 2 ^ 3 ). We therefore define the variable 
h by 

n - (Y - Y}a 2/ \ 

and in order to obtain some information about non-neutral disturbances, it is 
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convenient to work in a a“ 2 / 3 neighbourhood of 5^. It is of course 
possible to describe non-neutral modes for |x - XjJ » a -2 / 3 but it turns 
out that the Y variation of the disturbance is then on more than one scale 
so that the appropriate expansion procedure becomes much more tedious, for 
that reason we shall here restrict our attention only to |X - X | ~ a -2 ^ 3 . 
Thus, we define 

X = (X - X )a 2/3 
v n ' 


so that near (X n , Y n ) we can write 


< (X) = <Q + X< x a~ 2 ^ 3 + X 2 < 2 a 4/ ^ 3 + • •• 


u u Q0 + Xa U 1Q + na 2 / 3 u 01 + X 2 a + Xna 4 ^ 3 u n + n 2 a 4 ^ 3 u rt0 + 


02 


together with a similar expansion for w. The time scale for the instability 
now becomes 0(a 3 ^ 3 ), so a is expanded as 


8/3 J 2 4/3 

o = Cq a + Oj a + <*2 a ' + 


whilst the Gortler number G now expands as 


r _ „ 4 . 10/3 J 

G - g Q a + gj a + ••• . 


Finally, we expand the X velocity component U 


as 
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u = { u Q (x,n ) + a~ 2/3 u^x.n) + ••• } 


exp{ iot + a 2 / & Q (X) + a" 2/3 B^X) + ••• dx} 


together with similar expansions for V/a^, W/a^^, P/a®^. 

It now remains for us to substitute the above equations in the disturbance 

— 2/3 

equations (2.2) and to equate like powers of the small parameter a . 

However, before doing so, it is convenient to note that the operator 

* 8/3 

8/3t + u(3/3X) + iX a ' w expands as 


|_+ u |_ + X a 8/3 iw = a 4/3 {io 2 + u 0Q 3 2 + X + u Q1 n) 


+ 8 (/ U 20 X + u 02 ^ +u u Xh) 


+ iX (w 2Q X 2 + w Q2 n 2 + w x x 3bi )} + • * • 


a 4 ^ 3 M + ••• . 


Here we have assumed that B 1 are purely imaginary and satisfy 


ia 0 + P 0 U 00 + ^00 = °» 


ia l + 8 l u 00 + ^ W 10 X + 8 0 U 10 x °* 


6 0 U 01 + iXw 01 


0 


(2.14c) 


17 


The equations (2.14a,c) determine $q and ° q> whilst (2.14b) determines 
The second-order term in the expansion of o, namely c^, remains undetermined 
at this order; this is also the case at next order. Thus, o^, o 2 can 
apparently be chosen so that there is a continuum of three-dimensional modes 
having the asymptotic structure described above. 

It is now a routine matter to substitute the above expansions in (2.2) and 
equate like powers of a~ 2 ^. The zeroth-order system gives 

v o = g o K o u oo V 


together with the eigenrelation 


g 0 K 0 u 00 U 01 1# 


At next order we obtain a pair of linear equations for (Uj, V L ) which have a 
consistent solution if 


a 2 u. 


3n 


— + [o 0 + iY 0 ]U 0 + [Oj + iyJXUq + [a 2 + iY 2 ]nU 0 


+ iAy 3 x 2 U Q + iXy 4 xn U 0 + iXY 5 n 2 U Q = 0, 


(2.15c) 


where the real constants a Q , ctj, etc. are defined by 


°0 


U 00 g l K 0 2 

3 " 3 u 00 P 2r’ 


(2.16a) 


Y 0 _ + u 00 P 2i^ 


(2.16b) 
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a l 3 g 0 K ol- U 10 U 01 + U 11 U 00-l + 3 g 0 K l U 00 ’ 


2 g l u 10 

Yl " 3 u oo 


a 2 3 g 0 K 0^ u 10 + 2u 02 u 00-1 ’ 


2 U 01 


Y 2 = 3 °1 u 


00 


Y 4 = 


_ U 10 

W 01y 

w 01i 

0 

i — < 

1 O 
o 

1 3 

1 

" “o7 ' 

U 20 u Q1 J ’ 

U 11 

“01 , 

u 10 W 01^ 

' “01 “ 01 ' 


“ — -JJ* 

01 


Y 5 = 


2 r U 02 i 

3^02 ~ u Q1 W 0lJ * 


(2.16c) 
(2. 16d) 
(2. 16e) 
(2.16f) 

(2. 16g) 

(2. 16h) 
(2. 16i) 


Thus the vertical structure of the disturbance is determined by (2.15) and the 
most dangerous eigensolution is 


U 


0 


exp 




“2 + iy 2 Y 4 *** 
2(iXy 5 ) 2(iXy 5 ) 


2 

} 


whilst the corresponding eigenrelation becomes 

~ {a 9 + iY, + i^Y 4 x} 2 

{ ot Q + iY 0 } + {a 1 + iY^X + iXy 3 X Z ^ /-1 Xy 5 = 0. 


We have assumed that the disturbances are neutral at X - 0 so that by taking 
& 2 V = X = 0 the neutral value of gj 


is determined as 
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u 00 K 0 g l 

3 



a 2 Y 2 
2Xy 5 * 


Note that we have not needed to use the condition that u has a 
maximum at Y = Y in the above calculation. However, since we are interested 
in the most dangerous modes we now impose this condition in order to produce a 
minimum value for gQ. In this case = 0 and the above expression for 
gl simplifies to 


g i 




-I 

0 


in which case the flow is locally neutrally stable at X n if G expands as 


K(X n )u 00 U 01 ^ 


+ 3a 10/3 / 


x|y, 


U 01 + 


(2.17) 


Thus, the first correction term from the two-dimensional result is now 

0(a^^) and always has a stabilizing effect. If we replace X by in 

a 

(2.17), we obtain an asymptotic expansion of G in terms of X and a which 
remains valid until the two terms shown in (2.17) are comparable. This occurs 
when X is formally O(a^) and this is the regime we now consider. We shall 
show that at this stage the neutral value of the Gortler number can be shifted 
by an 0(1) amount from its two-dimensional value. Moreover, we shall see 
that at this stage the eigenfunctions essentially reduce to Orr-Sommerf eld 
"centre modes" destabilized by centrifugal effects. 
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_ 3 4 

c. The limit a > ~ with X ~ 0(a ), G ~ 0(a ) . 

We first note that if we wish to retain the critical layer structure of 

— A *3 / 

(b), with X now replaced by X a , we must formally take 3/3X~ 0(aj in 
(2.2). It then follows from the continuity and momentum equations that this 
is possible only if W ~ a 3 , V ~ a 2 , P~ a 3 , with U again taken to be 

0(1). Let us again assume that the location of the critical layer at the 
neutral value of X = is Y = Y and define stretched variables n 
and X by 

n = {Y - Y)a, X = {X - X n >a, 

so that the vortices have their depth and wavelength comparable. The 
disturbance then expands as 

A,/ 

u = { Uq (n ,X) + a 1 l^Oi.X) + •••} exp{iat + a 3 / $ Q (X) + a -1 P^X) + ••• dx} 

together with similar expansions for V/a 2 , W/a 3 , and P/a 3 . We retain the 
expansions of k, u, w given in (b) but with a“ 2 / 3 replaced by a -1 . 

Finally, the frequency and Gortler number expand as 

4 3 

a = a Og+a o ^ + • • • , 

^ 4 3 

G = a g Q + a gj + • • • . 

The eigenrelations (2.14a,b,c) are again taken to be satisfied so that 

A / 

3/3t + u(3/3X) + Xia w expands as a 2 M with M defined by (2.13). 
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We first consider the equation of continuity, which at order a\ a^ 


yields 


*0 u o + iw o = °’ 


9v o 

^o u i + iw i = - e i u o"— • 


The second of these equations can be written in the form 


3V ( 

iW l =“ e i U 0 - W 


(2.18a) 


if we first define = ip^ + W^. The three momentum equations then 

yield 


LU o = U 01 V 


(2.18b) 


3P 0 

LV 0 3n + g 0 K 0 u 00 V 


(2.18c) 


LW X = iP Q + XV Q {2y 6 n + Y i X} , (2. 18d) 

where T 6 = 2{» 02 ‘ ^ V • r 7 " W 11 ‘ ^ »01- and with 

M defined by (2.13). The system (2.18), together with the conditions Uq, 

Vq» W - 0 when |n | -*■ 00 , specifies an eigenvalue problem for g + i0 2± 

■* ~ 

if X, a 2 » and X are given. For simplicity, we consider only the neutral 

case g„ = 0 so that we can set X = g„ =0 in (2.18). Furthermore, we 

2r 2r 

will discuss the solution of (2.18) only for a \ ~ a 2 = we not P rove 

A 

that this will lead to the most dangerous modes for all X , but it is 
certainly the case for X « 1. The system (2.18) can then be simplified to 
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{ 2 “ 1 “ u 00 ^2 “ X * Y 6 n 2 " ^ V 0 + 2l * Y 6 V 0 ” g 0 K 0 u 00 U 0 ’ ^ 2 * 19a ^ 



U 00 3 2 



*01 V 


(2.19b) 


Here $ 2r = 0 and 3 2i , must be chosen such that (2.19a,b) are satisfied, 
and with Uq, Vq + 0, n -*• if°. It is clear from (2.19) that gp(-A) = gQ(A), 

A A 

0 2 (-X) = P 2 (X) so without any loss of generality we shall from now on assume 
that Ay^ > 0. It is convenient at this stage to note that the eigenrelation 
for gg can be written 

g 0 K 0 u 00 U 01 = -^ Y 6 ( 2 . 19 c ) 


where A) must be calculated numerically. However, if we let A -*■ 0 

in (2.19), we obtain 

* A 2/3 

^(Y 6 A) ■* 1 + 0(A) Z/J , 


which is consistent with the results of (b). 


Thus for a » 0(A 1 ^) 


g 0 K 0 u 00 U 01 ~ ** 

which is, of course, the two-dimensional result* When y^ X ■* » , it can be 

shown that the functions Uq and Vq have a double boundary layer structure 

* -1 /4 

corresponding to n = 0(1) and n = O(Ay^) • In this limit we write 

c(Ay 6 ) 3/2 + ••• 

u 00 P 2 = id < AY 6 )1/2 + 
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a j [t 

and define £ = (Ay^) n . In the inner region where t, = 0(1), we write 


u 0 =Uo<0 + °(iY 6 r 1/2 , 


~ «V‘ /2 ~o 

V = TT 

0 u oi 


Vq(?) + 0(1), 




and then U Q , V Q satisfy 


{!_ - id - i? z } 


9? : 

j 9 ^ 2i~0 _ ~0 

{ 2 d * U 0 " V 


9 ~ 0 ~ 0 ~ 0 

1 y + 2i V = -c U u 

v 0 + Z1V 0 C V 


( 2 . 20 ) 


In the outer region, where n ~ 0(1), the first approximation to (2.19) 
reduces to a second-order differential equation for Vq independent of Uq 
and a linear equation which involves Uq and Vq. The solution of these 
equations in terms of Bessel functions shows that the condition Uq, Vq 0 
when X] °° leads to the following matching conditions on (2.20): 


V 


0 

0 


/V 


5 ’ 



c * °° . 


(2.21) 


The most dangerous modes of (2.20) are even in £ so we solve (2.20) subject 
to (2.21) and 

/■ ^ ^ ^ ^ 

0 

Uq = Uq = Vq = 0 iC - 0. (2.22) 


The eigenvalues c,d associated with (2.20), (2.21), and (2.22) were found by 
a shooting procedure using (2.21) and the two other independent but 
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exponentially decaying solutions of (2.20) which can be found for |? | >> 1. 
The three independent solutions were integrated from a suitably large value of 
C to x, = 0 where they were combined to satisfy (2.22). The least stable 
mode corresponds to 

c = 4.71, d = -2.89. 


For intermediate values of and 0 2i were calculated numerically by 

solving (2.19) using a shooting procedure with the three independent exponen- 
tially decaying solutions of these equations. Figures 1 and 2 illustrates the 
results of such a calculation for the least stable mode which is an even func- 
tion of ii. In these figures we have also shown the one-term asymptotic 

A 

approximations to the neutral curves valid for y^ \ » 1. The monotonic 
increase of the function yV has a profound and unexpected influence on the 
nature of the neutral curve. We shall now show that for X » 1 a neutral 

curve with a minimum value for the Gortler number is predicted by our 
calculations. 

In order to see why this is the case we first note that (2.19c) written in 
terras of G and X gives 

<*o u oo U 01 = » (2,23) 

a 

and the previously described asymptotic structure for implies that 


<*0 u 00 U 01 ~ a ’ 


X 1/3 « a, 


4.71(Xy 6 ) 3/2 /3 

<*o u oo u oi T72 » x >> a - 
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Thus if A is held fixed, the Gortler number decreases or increases with the 

—1/3 

wavenumber, depending on the size of a/X . Hence, there will be a minimum 
value of G at some 0(1) value of a/A^^. In fact, it is easily shown 
from (2.23) that this occurs when 

^'(Y 6 X) " — X), 

3y 6 X 6 

A 

and our calculations show that this occurs where y, X = 3.18 and the 

b 

minimum value of G is given by 

G K o u oo u oi = 3,88 (y 6 ^ » 

and then 

a = .68 (y 6 

In Figure 3 we have shown the dependence of G k q u qo u Q1 on a implied 

by (2.23) for several values of Xy^. We note that for a particular three- 

dimensional boundary layer the quantities k q> u qo , u q1 ,y 6 must be calculated 

for a particular choice of (X n ,Y) and then (2.23) and Figure 1 can be used 

to generate the neutral curve for the required values of T, The crucial 

result shown in Figure 3 is that for |T| » 1 a minimum Gortler number 

exists for a disturbance concentrated at a fixed height above the wall. The 

eigenfunctions associated with some of the results of Figures 1 and 2 are 

shown in Figures 4 and 5. We note that the vortices spread further away from 
/\ 

the wall when Xy^ increases. Suppose then that we have a three-dimensional 
boundary layer flow with X » 1 but fixed. For any given height T of the 
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—4/3 

critical layer we have shown that neutral disturbances occur for G > 0(A ) 

with a ~ This means that the Gortler mechanism in a three-dimensional 

boundary layer with X >> 1 is not operational at 0(1) values of the 
Gortler number. Since the crossflow on a laminar flow wing is certainly large 
compared to R -1 / 2 , with R scaled on the chordwise velocity, we have 
effectively shown that the Gortler mechanism is almost certainly unimportant 
in such a flow. Thus, it appears likely that crossflow instabilities 
associated with inflection points in the basic flow velocity component in 
particular directions will be most important in determining the state of the 
boundary layer. 

It is interesting to note that if Uq is set equal to zero in (2.19), 
then the equation for V Q reduces to the rescaled Orr-Sommerfeld equation for 
a two-dimensional Tollmien— Schlichting wave, having its critical layer located 
where the velocity field has a local maximum. Such 'centre modes' are, of 
course, stable and have a wave-packet structure in the Y direction similar 
to that discussed by Tatsurai and Gotoh (1969). Thus, we interpret the neutral 
modes found above as Tollmien— Schlichting waves destabilized by centrifugal 
effects. 


3. THE DEGENERATE CASE OF FLOWS WITH ZERO PRESSURE GRADIENT 

The analysis of Section 2 does not apply if u and w are linearly 
related as is the case when the pressure gradient in the boundary layer 
vanishes. Suppose that this is the case and the basic flow is written in the 


form 
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u = (u(X,Y),R _1/2 v(X,Y) ,cu(X,Y)) ( 1 + 0(R R)) , 

where c is a constant and X,Y are as defined in Section 2, whilst u and 

V are determined in terms of the Blasius function f( . We perturb the 

/X 

above flow such that the disturbance velocity field is periodic in the X 
and Z directions, but with the wavefronts parallel to the basic flow. Thus 
we write the perturbation velocity field u' as 

= (UE, R -1 ^ 2 VE, cUE + R _1 ^ 2 WE(1 + c 2 ) 1 ^ 2 )(l + 0(R _1 ^ 2 )) 

U 0 

where E = exp i[aZ - Xa] . The 0(1) contribution to the Z component 

of the disturbance velocity field has been inserted in order that the 

component of the disturbance velocity field in a direction perpendicular to 

the basic flow direction is 0 (R~ iy ^). we also assume that the perturbation 

is independent of time and the pressure perturbation P(X,Y) corresponding to 

-1 2 

the disturbed flow is scaled on R p Uq. In order to determine the 

equations satisfied by U, V, W and P, it is convenient to use the X,Y 

momentum equations and the Z momentum equation subtracted from the X 

momentum equation multiplied by c . The particular X-dependence of the 

perturbation we have chosen means that 9/9X in these equations and the 

1/2 9 

continuity equation should be replaced by -cR ia + ^ • After some 
manipulation it can be shown that the equations satisfied by U, V, W and P 

U x + V Y + iaW /l + C 2 = 0 

lST x + uU x + vV y + Vu y = U YY - a 2 U[l + c 2 ] 


are 
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uV X + Uv X + vV Y + Vv Y + Gk(X > uU = “ P y + - a 2 V[1 + c 2 ] 

^W x + ^y “ ~ iap f 1 + c 2 ] 1/2 + W yY - a 2 W[1 + c 2 ]. 

The above equations are in fact identical to those appropriate to the two- 
dimensional boundary layer flow (u, R *^ 2 v, 0) if we identify Jl + c 2 a 
as the 'total' wavenumber of the perturbation. Hence, the two-dimensional 
results of I, II are applicable and no further calculation is necessary. We 
note, however, that in this case, the vortices are aligned with the basic flow 
and that the effective Gortler number is determined by the component of the 
basic flow normal to the leading edge. 


4. CONCLUSION 

Let us first consider the implications of the results of the last section 
which are appropriate to boundary layers with zero pressure gradient. Suppose 
then that we have a Blasius boundary layer over a wall of variable curvature 
with the leading edge perpendicular to the flow direction at infinity. 
The Gortler vortex instability mechanism described in I, II is operational in 
such a flow and the manner in which the instability develops in the streamwise 
direction depends on the Gortler number G and the wavenumber a. If the 
wavenumber a is not large, this development also has a strong dependence on 
the form of the initial disturbance and the concept of a unique neutral curve 
is not tenable. Suppose next that the leading edge is inclined at an angle 
a to the oncoming flow. The basic flow again consists of a Blasius boundary 
layer with the velocity field in the plane of the wall parallel to the flow 
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direction at infinity* We choose to look for a Gortler vortex instability 
with the vortices parallel to the flow at infinity and define the Gortler 
number with respect to the component of the flow velocity at infinity 
perpendicular to the leading edge* The disturbance equations which we found 
in Section 3 are identical to those appropriate to the two-dimensional flow 
obtained by neglecting the crossflow. Thus, when the leading edge is turned 
through an angle a, the vortices remain parallel to their original direction. 
However, their development in the strearawise direction will vary with a, 
since the instability is governed by the Gortler number scaled on the velocity 
component perpendicular to the leading edge. 

Let us now turn to the results of Section 3 which correspond to the flow 

- 1/2 

over a wall swept at angle a ~ 0(R ) to the oncoming flow. The results 

we have found are quite general and it was not necessary for us to be specific 
about the form of the basic boundary layer flow. The results apply to, for 
example, the situation where the zero sweep solution is given by the Falkner- 
Skan solutions or to the case when the basic three-dimensional flow is given 
by the Sowerby-Loos solution (see Rosehead (1963)) which is, of course, 
relevant to the flow over a rotating blade. 

The results of Section 2 were obtained by formally taking the crossflow 
velocity field to be 0 (R“*^) and then considering various asymptotic limits 
involving a, G and X. We have concentrated our attention on the limit 
X + °°, since we are interested in finding the structure of the instability for 
crossflow velocity field large compared to R” ^ ^ . We have shown that for a 
three-dimensional boundary layer, the most dangerous modes are time-dependent 
and the orientation of the vortices in the neutral location is fixed by the 
vortex lines of the basic flow and not by the streamlines, as we would expect 
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based on the results for the two-dimensional problem. The most surprising 
result we have found is that for the three-dimensional problem with X » 1, 
the large wavenumber asymptotics developed in I produces a neutral curve with 
distinct left- and right-hand branches and a minimum Gortler number. Each such 
neutral curve corresponds to a disturbance with its critical layer fixed at a 
given height above the wall. 

In essence we have shown that, since the minimum value of the Gortler 
number on any such curve is 0(X"^^), the Gortler vortex instability mechanism 
is almost certainly unimportant in the concave region of laminar flow swept 
wings where the angle of sweep is large compared to R • Therefore it 
would appear that crossflow or Tollmien-Schlichting wave instabilities (which 
occur for Reynolds number typically of order 100) are more likely to lead to 
the breakdown of laminar flow. 

We can, in fact, be more precise about the angle of sweep at which the 
Gortler vortex mechanism ceases to be important. This can be done by noting 
that the Gortler number was initially defined in terms of the two independent 
parameters 6 , and R, so that G >)> 1 implies that 


6 » R 


Thus, the analysis of Section 2c can be repeated by writing the disturbance 
equations in terms of R, 6 and X but with 

6 ~ R _P , p > 0 

X ~ R +q , q > 0. 
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The whole expansion procedure can then be reformulated in terms of R as the 
large parameter, and we can determine for what range of values of p and q 
the exact structure of Section 2c is recovered. Such a calculation shows that 
the results of Section 2c remain valid for 

q < 3/8, 

so that the critical Gortler number is asymptotically large whenever the 
crossflow is 0(R There seems little doubt, therefore, that for 0(1) 
angles of sweep, the Gortler vortex mechanism is not significant in the 
Reynolds number regimes where crossflow and Tollmien-Schlichting instabilities 
occur. 
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Figure 1* The dependence of jV on Ay^, the dotted curve corresponds to the 
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one-term asymptotic expansion of ^ valid for Ay. » 1. 
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